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The quantum anomalous Hall system with Chem number 2 can be destroyed by sufficiently strong disorder. 
During its process towards localization, it was found that the electronic states will be directly localized to an 
Anderson insulator (with Chern number 0), without an intermediate Hall plateau with Chern number 1. Here 
we investigate the topological origin of this phenomenon, by calculating the band structures and Chern numbers 
for disordered supercells. We find that on the route towards localization, there exists a hidden state with Chem 
number 1, but it is too short and too fluctuating to be practically observable. This intermediate state cannot be 
stabilized even after some “smart design” of the model and this should be a universal phenomena for insulators 
with high Chern numbers. By performing numerical scaling of conductances, we also plot the renormalization 
group flows for this transition, with Chem number 1 state as an unstable fixed point. This is distinct from known 
results, and can be tested by experiments and further theoretical analysis. 

PACS numbers: 


Introduction 

As well as integer quantum Hall effect (QHE) under a magnetic fieldi, the quantum anomalous Hall effect (QAHE) without 
an external magnetic field^ is characterized by nonzero Chern number, a topological invariant associated with occupied bands^. 
This nonzero Chern number C gives rise to edge states in the bulk gap, which carry the remarkably perfect quantization of 
the Hall conductance cr^y - C (in units of ^ hereafter). Inspired by the investigations of topological insulators^2, QAHEs in 
concrete materials have been theoretically proposed^^ and been experimentally observecU^— . 

Although robust against weak disorder, topological phases will be localized into Anderson insulators ultimately, by sufficiently 
strong disorderi^— . Microscopically speaking, this can be viewed as a disorder induced topological phase transition, i.e., change 
of the topological invariant (1 —> 0) due to a band touching, i.e., gap closing and re-openin g^*'^^ . 

Topological phases with different Chern numbers are topologically distinct from each other. Each of them corresponds to a 
stable fixed point in the sense of renormalization group (RG)^, which offers another explanation for their robustness. Numerical 
calculations for a QHE with a high Chern number C > 1 show that, with the increasing of disorder strength, the Hall conductance 
cr^y - C vanishes to 0 persistently without showing any intermediate Hall plateaus with Chern numbers C - 1, C - 2, • • • , 1—^, 
which are also associated with stable fixed points and therefore should be robust. The quantum Hall state with C > 1 is ascribed 
to contributions from individual occupied Landau levels with C = 1. Disorder makes the Landau levels broaden. After adjacent 
Landau levels are broadened enough to touch each other, their Chern numbers change from ±1 to 0, which results in trivial 
localized bands. The most essential observation is that, this disorder induced topological transitions happen to all Landau levels 
almost at similar disorder strength, which means a nearly simultaneous breakdown of all Landau level a*^'^* . Therefore, after 
disorder ensemble average, any intermediate topological state with only a fraction of the Landau levels surviving cannot stably 
manifests itself. 

Lattice models for QAHE with a single valence band carrying a high Chern number C > 1 were proposed recently2i^. 
The localization process with increasing disorder strength for a two-band model with C = 2 was investigated by transport 
calculations2i. It was found that its Hall conductance still decays from 2 persistently to 0 without a plateau 1. Now since one 
single band carries C = 2, the above picture of simultaneous breakdowns for all Landau levels (each carrying C = 1) for QHE^ 
does not apply. Without special symmetry protection, the band touching induced by a single parameter (disorder strength here) is 
linearly shaped and only changes the Chern number by Therefore, the topological origin of above phenomenon seems to 

be confusing. Does this correspond to a topological phase transition from two band-touching points happening simultaneously, 
or from a single band-touching point changing the Chern number by 2? 

In this manuscript, based on the language of disordered supercell s^'^* , 22 , 30,3 1 ^ investigate the topological origin of the 
localization process for QAHE with Chern number 2 —> 0. We find that, with the increasing of disorder strength W, this transition 
corresponds to two successive band-touchings which happen at different disorder strengths Wi and W 2 , and at different positions 
fci and k 2 in the Brillouin zone (BZ). The intermediate C = 1 window (Wi, W 2 ) is small and highly configuration dependent. 
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FIG. 1: (a) The band structure of the ribbon in clean limit with width L,. = 100. (b) Two-terminal conductance as a function of disorder 
strength, which is averaged over 500 disorder configurations with sample size = Ly = 100. The model parameters as follows: t\ = 1.0, 
t 2 = 1.0, f 3 = 0.5, OT = -1.0, and Ef = Q.l. 



which makes it hardly observable after disorder averaging and size scaling. This vulnerability of intermediate C - 1 state cannot 
be improved even if we generalize the original model to make the gap at one valley much smaller than at other valleys (thus 
with the expectation that the band-touching at the smallest gap should be much earlier). Finally, based on finite-size scaling of 
mesoscopic conductances, we summarize the results into a diagram of RG flows, which can be tested by experiment a^^’^^ or held 
theoretical analysis^^. 


Results and Discussions 

We use a generalized version of the two-band model recently propose d^^’^^ . This is a spinless model defined on a square 
lattice with two orbitals at each site. Setting the lattice constant a = 1, the Hamiltonian in k space reads^i 


Hoik) - dia,c +d2cry + ( 1 ) 

d\{k) - 2t\ coskx, d 2 {k) = 2fi cosky, 
d^ik) - m + 2tj, sin kx + 21^ sin ky + 2 t 2 cos{kx + ky), 

where cr, are Pauli matrices in the orbital space. After realizing the Hamiltonian on the square lattice, one can see that m is the 
staggered potential on two orbitals, and f, (i + 2) are the nearest hoppings between different sites. The next nearest hopping 
f 2 just exists along one direction, e.g. along two-dimensional crystallographic orientation [11], which is essential for realizing 
Chern number 2r^^. Hereafter, fi = 1 will be used as the energy unit. In the BZ, model ([TJ has four Dirac points (or valleys) at 

D\XT>A - {-j’ ±f) associated masses (or gaps) A,- are 2((m - 2t2) ± 4f3) and 2(m + 2t2), respectively. 

The Hall conductance is related to the totol Chern number of all occupied bands as cr„ = x Chern number 

associated with the n-th band is defined in two-dimensional k-space: 

r d^kF,{k\ (2) 

Am Jsz 

F = X A, A = <n(fc)|V,t|n(fc)), 

where F is the Berry curvature, and |n(A;)) is the normalized wave function of the n-th Bloch band such that H{k) \n(k)} = 
En{k) \n(k)}. For our model with Hamiltonian ([T]), the Chern number associated with the valance band is straightforwardly 
determined by the model parameters as^S. 

C - -[sgn(m - 2t2 + A-t^) + sgn(m - 2t2 - Tfa) - sgn(m + 2t2) — sgn(m -i- 2f2)]. 

The effect of disorder is included in the real space representation by adding a random potential term t/,- to each site i, with t/, 
uniformly distributed in the interval {-WI2, W/2), where W is the disorder strength. The Hall conductance of a disordered sample 
can be simulated by using the standard method of non-equilibrium Green’s functions^. On the other hand, disorder breaks the 
translation symmetry of the original lattice and makes the original k space badly defined for Chern number in equation (|2]i. 
Nevertheless, for a disordered sample with size Ly x Ly, if twisted boundary conditions exp{ikxLx) and exp{ikyLy) are introduced 
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FIG. 2: (a) Three quantities as functions of disorder strength W for a definite disorder configuration: the hand gap A (red line), the Berry 
curvature F (green line) and the Chem number C of the valence band (blue bar), (b) The location of minimal gap point (red circles) as a 
function of disorder strength W for the same disorder configuration. The cyan lines indicate their projections on the k,, - ky plane, (c) The band 
gap A (red line) and the Chem number C below the gap (blue bar) as functions of disorder strength for 8 different disorder configurations. The 
intermediate windows (ITi, IV 2 ) with C = I are marked by green arrows. 


in both directions, the definitions of fe = (ky, ky) and Chern number can be restoredi^. This is not surprising since this is 
equivalent to constructing a superlattice with this Ly x Ly sample as the supercell22i2i. Physically reliable results for “really” 
disordered systems will be recovered after disorder ensemble averaging and size scaling L —» 00. See Methods for details of 
these calculations. 

We first reproduce the main results of previous transport calculations^i. Fig. [1] (a) is the band structure for a quasi-one 
dimensional ribbon, where two pairs of edge states from C = 2 can be seen. Although there are essentially 4 massive Dirac 
points for the two-dimensional band structure, only 2 of them are distinguishable here due to the cutting orientation of producing 
the quasi-one dimenional ribbon. In Fig. [I](b), we display the the averaged two-terminal conductance G (- cr^y when C 4 0) as 
a function of disorder strength W . The conductance of the C — 2 state is perfectly quantized regardless of finite disorder, until 
W ~ 6. Afterwards, as found previously, the conductance decays from 2 towards 0 persistently, without any plateau at the C = 1 
state (the dashed line)S. 

In order to reveal the topological origin of the above transport results, we resort to the calculations of band structures and 
Chern numbers by using the disordered supercell method introduced in Methods. Random potential Ui - W ■ m, are assigned to 
each site i of a supercell with L^xLy sites, where m, are random numbers uniformly distributed within the unit interval (-0.5,0.5), 
and W is the disorder strength. Hereafter, {f/,} - W ■ Ui from a fixed configuration of {m,} but only differed by the factor W will 
still be called one disorder configuration (or one sample). Applying twisted boundary conditions to this disordered supercell, we 
can investigate the adiabatic changes of band structures and Chern numbers with the increasing of the single parameter W, for a 
fixed disorder configuration {m,}. 

Let’s first focus on the microscopic evolution of one typical disorder configuration. To find possible band touching, we define 
the band gap A as the minimal difference between conduction and valence bands £’c(v)(fc) (with the same k) in the BZ, 

A s minjFcCfc) - Ey(k)}, 

keBZ 

with the associated point kyy^in as the minimal gap point. A band touching corresponds to A = 0. At the clean limit, kmin is just 
one of the Dirac points D, . In Fig. |2](a), we plot the band gap A (red line), the Berry curvature F at the minimal gap point kyyyiyy 
(green line) and the Chern number C of the valance band (blue bar) as functions of disorder strength W, around the transition 
region. It is interesting to notice that for a definite configuration, there are two subsequent topological transitions at Wi ~ 7.2 
and IT 2 ~ 7.7, with band touchings (A = 0), and with Chern numbers 2 —» 1, then 1^0, respectively. In Fig. |2](b), we plot the 
track of the minimal gap point kmin in the BZ as a function of disorder strength W. Approaching the first topological transition 
Wi ~ 7.2 from below, begins to digress from Dirac point Di = (|, |). Soon after this, kmin appears around another Dirac 
point D 3 = (-J, -j), and the second topological transition Wi ~ 7.7 happens around here. Such large transfer of the minimal 
gap point is not surprising since the topological phase collapses at strong disorder, which has remarkably changed the original 
valley shapes in the clean limit. 
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FIG. 3; The same with Fig. [T]but with model parameters as ?i = 1.0, t 2 = 1.0, ts = -0.72, m = -1.0, and ii/r = 0.1. (a) The band structure of 
the ribbon in clean limit, (b) The disorder averaged two-terminal conductance as a function of disorder strength. 


After confirming the existence of an intermediate C = 1 state for a definite disorder configuration, we have to answer why it 
cannot be observed in the transport calculations in Fig. [T|(b). In Fig. |2](c), we show the gaps A (red curve) and Chern numbers 
C (blue bar) as functions of W for 8 different disorder configurations, each with an intermediate C - 1 state between Wi and IF 2 
(the green arrows). The first observation is that these C = 1 windows are rather narrow, i.e., W 2 -W 1 <& Wi (Note the leftmost of 
the VF-axis is 6 instead of 0.). Another important observation is that such windows {W\, W 2 ) with C = 1 are quite configuration 
dependent, with their locations and lengths fluctuating remarkably from configuration to configuration. There is not a stable 
common window with Chern number C = 1 for all different configurations. Furthermore, we have checked that, this absence 
of a stable common window for C = 1 does not get better at all, with the increasing of sample size L. In short, although there 
exists an intermediate window of C = 1 for each definite configuration, this window is too vulnerable to survive after disorder 
averaging and size scaling, and thus is physically unobservable. This is the first important finding in this manuscript. 

The two subsequent band touchings at Wi and W 2 are too close to make up a stable state with C = 1 within (Wi, VF?)- One 
may blame this to the almost identical gap sizes at two Dirac points in the clean limit of the model (See Fig. [T|(a)). This tempts 
us to modify the model to make the gap at one Dirac point, say, Ai much smaller than those at other Dirac points. This can be 
realized for model ([1]) by setting parameters as depicted in Fig. [3] Then, under perturbation from rather weaker disorder—, we 
expect the first band touching VFi to occur around this Dirac point Aj with the narrowest gap, much earlier than W 2 occurring 
around another Dirac point. Thus a wider window (Wi, W 2 ) could be expected and this could be beneficial for stabilizing the 
C = 1 even after disorder averaging. 

In Fig. [3] (a), the band structure for a ribbon is plotted. The gap at Di is 0.24, only 4% of that at £> 2 . The two-terminal 
conductance as a function of increasing disorder is plotted in Fig[3(b). Compared to Fig. [I](b), the C = 2 plateau starts to 
collapse earlier at IF ~ 4, which can be attributed to the small gap at the clean limit. However, a stable plateau of C = 1 is still 
absent. 

We repeat the calculations of disordered supercells to reveal the underlying topological nature. In order to see the evolutions 



FIG. 4: With model parameters same as in Fig. [3] the developments of quantities for a disorder configuration as functions of disorder strength 
IF: the local gaps at 4 Dirac points A, (solid lines), the minimal band gap A (red dotted line) and the Chern number C (blue bar). Note that 
after W > 3.6, the minimal band gap will not be located at any of the Dirac points any more. 
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FIG. 5: Distributions of (crjcy for 10000 disorder configurations with supercell size 20 x 20, for different disorder strengths W. The Fermi 
energy is fixed at Ep = -0.053, near the gap centre of the clean limit. Other model parameters are same with Figs.[T] 


of different massive Dirac points, local gaps at 4 Dirac points D, are simply defined as 

A, = E,{Di) - 

Without disorder, Ai is just the band gap A, which is defined at the narrowest point in the gap. In Fig. 01 we plot the developments 
of A, (solid lines), as well as A (red dotted line) when increasing disorder strength W for a typical disorder configuration. The 
Chern number (blue bar) is also plotted as a reference. In clean limit W = 0, as we set, Ai (green line) [which is equal to A 
(red dotted line)] is much smaller than other 3 local gaps. With the increasing of W, those 3 large gaps decrease. Contrary to 
what we are expecting, the smallest gap Aj is increasing to prevent itself from an early band touching at weak disorder. After 
W ~ 3.6, A does not follow Ai (nor any other A,) any more, which means the minimal gap point has drifted away from the 
Dirac points. This, as discussed above, is a consequence of remarkable destroy of the original shape of the band structure from 
strong disorder. In other words, the elaborately designed band structure with extremely small local gap at one Dirac point Ai 
has been smeared out completely by disorder, before the first band touching occurs. Then, as happened in Fig. |2] there will be 
two subsequent and close band touchings, which fluctuate so seriously that the C - 1 state between them is unobservable after 
disorder averaging and size scaling. In one word, in the BZ, states around all 4 Dirac points seem to “feel” each other and tend 
to evolve cooperatively under disorder. This reflects the fact that the Chern number C and its transition are global properties of 
the band^— , especially in the presence of disorder, i.e., without constraints from most spacial symmetries. We conjecture that 
for Chern insulators with high Chern numbers C > 1, either QHEs with multiple Landau levels or QAHEs with a single valence 
band, this route towards localization under disorder with Chern number C ^ 0 is universal, independent of the concrete details 
of the materials in the clean limit. 

It is interesting to revisit this route towards localization, by observing the evolution of the distribution of the Hall conductance. 
In Eig. |5] we present the statistical histogram of cr^y at a definite Eermi energy and sample size, for different disorder strength 
W. We will see that they also reflects the physical pictures described above. For W - 6 when (Txy starts to decrease from 2, Fig. 
|5](a) shows that a small portion (around 10%) of samples have experienced a topological transition 2 —> 1. For a larger W as 
illustrated in Fig. |5](b), besides two peaks at cr^y - 1,2, a considerable amount of samples with non-quantized (Txy appear. They 
originate from partly filled subbands randomly distributed around the Fermi energy, which are induced by strong disorder. Such 
a broad distribution gives rise to a non-quantized cTxy with a large statistical variation. When W = 10 [Fig. |5](c)], most C - 2 
states have collapsed. However, there is a prominent and sharp peak at C = 1. This means that Ef - -0.053 here is within the 
(small but finite) gaps of these samples with C = 1 below them; Or in other words, VT = 10 here is within the C = 1 windows 
{W]_, Wi) [e.g., those marked by green arrows in Fig|2](c)] of these samples. Nevertheless, as stated earlier, there is no common 
window with C = 1 for all (or most) samples. For the rest samples (more than 80%), their cr^ are still widely distributed, 
making a stable and observable cTxy - 1 state impossible. Fig|5](a) to (c) reaffirm the “hidden” nature of the intermediate C = 1 
state: it is existent but overwhelmed by strong statistical fluctuations, and thus hardly observable. This is contrary to the case of 
topological Anderson insulator—. In that case, although also with dense disorder induced in-gap states, the topological invariant 
still has a probability of almost 100% to be 1, in a finite range of energy and disorder strength, giving rise to an observable 
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FIG. 6: The renormalization group flows extracted from finite size scaling for square samples with = 10 —> 15 —> 20, with Fermi energies 
Ep 6 (-1,1) and disorder strengths as shown. Other model parameters are same with Fig.[T]and each point is an average over 10000 disorder 
configurations (50000 configurations for W = 9, 10). The dashed line is the semi-ellipse with the semi-major (semi-minor) axis a = 1 (b = cti), 
where <ti = |. The C = 1 state corresponds to point (1, cti). 


quantized conductance plateau with tiny fluctuations^. At the end of this transition, as shown in Figl5](d), most samples have 
crossed the second topological transition W 2 towards trivial localization with C -0. 

Finally let’s put this localization process into an RG flow diagram by numerical scaling. RG analysis can provide vivid and 
deep understanding of phase transitions. Besides, recent experimental progresses have made RG flows for QAHEs directly 
observablo^^i^i. The sample size driven (cr^^, cr^.^) flows are plotted in Fig. |6l corresponding to the localization process in Fig. 
[1] Similar to the RG diagram in the transition with C - 1 —> most flow arrows reside on a semi-ellipse. On two ends 

of the semi-ellipse, scaling flows are attracted to two fixed points (cTxy, (Txx) - (C, 0), corresponding to the C - 2 {C - 0) state 
before (after) the disorder induced transition, respectively. On the other hand, the intermediate state C = 1 is located on the 
uppermost of this semi-ellipse, with the maximum longitudinal conductance cti = 0.79 ~ 7 r/4. Around this point, some scaling 
flows are running away, and some are randomly oriented. The latter behavior can be attributed to the statistical precision limited 
by the number of disorder configurations we can achieve, but also reflects the unstable nature of this point. Thus it’s reasonable 
to conclude that this point (1, cri) is a saddle type unstable fixed point, like the role that (5, cr*) plays in the transition with Chern 
number 1 —> o 23,36,37 _ brief, we show that the state C = 1 does not always correspond to a stable fixed point in the sense 
of RG. Furthermore, the associated maximum longitudinal conductance cri = 0.79 ~ 7r/4, which is twice as that found in the 
localization process with Chern number C = 1 ^ 0 stated. These features make Fig. |6]distinct from known RG knowledge for 
Chern insulators so far. This can be tested by further theoretical analysis and experiments. 

In summary, by investigating the band structures and Chern numbers of disordered supercells, we illustrate the route towards 
localization for QAHE with C = 2. With the increase of disorder strength, there are two successive topological transitions 
for each disorder configuration, but the intermediate C = 1 is remarkably short and fluctuating, therefore cannot survive after 
disorder averaging and size scaling. This picture is also supported by statistical distributions and RG flows of the conductances. 


Methods 


At zero temperature, the two-terminal conductance of a finite sample can be expressed by Green’s functions ai 

G= -Tr[FsGTDG"], 
h 


,34 


( 3 ) 
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where G''^“ is the retarded/advanced Green’s function, and rs(/j) = ~ ^s(dO ^s{D) retarded/advanced self en¬ 

ergies due to the source (drain) lead, respectively. For a quantum (anomalous) Hall state (C + 0), this two-terminal conductance 
is just the Hall conductance, G - Ux^- The transport results in Figs. [I](b) and[3](b) were calculated by this method. However, 
for topologically trivial state, e.g., a state close to localization, this simple two-terminal G may not rigorously equal to (Txy since 
it will contain bulk transports. Another method of calculating cTxy will be introduced below. 

Numerical methods are necessary for treating disordered supercells with large size L^xLy. However, the numerical calculation 
of topological quantities is a tricky task: Special care must be taken to ensure the gauge invariance in the numerical processes. 
We adopt a method of discretizing equation (|2]) as^ 


c„ = ^^7-(k,). (4) 

Here / runs over all lattices over the discretized BZ. On this discretized BZ, the Berry curvature 'F{ki) is a lattice field defined 
on each small plaquette as 'F{ki) = \n'Ll\{ki)^ 2 {ki -i- l)'L(i(ki -i- 2 )“*'Z/ 2 (k;)“', where ^fi(ki) = 1 and 2 are unit 

vectors of of the discretized BZ^. 

For the purpose of plotting RG flows, both the longitudinal conductance cr^x and the Hall conductance cr^y of the disordered 
sample are needed. A simulation of six-terminal Hall bar by using non-equilibrium Green’s functions like equation will 
be numerically expensive. Instead, we adopt the methods which were successful for plotting RG flows^. The Hall conductance 

CTxy is calculated by a numerical integration of Berry curvature under the Fermi energy Ef using equation Q. For a given Ef, 

2 

(Txy may not be quantized if some subbands are partly filled. The longitudinal conductance cTxx (also in units of ^) is calculated 
as the Thouless conductance ^^4° 


x(E„) = 7Tp{E„) 


^1 

dki I <=1=0 


(5) 


with twisted (periodic) boundary condition in x (y) direction respectively, where p is the density of states, and £’„(k) is the n-th 
subband. The conductances in Figs. |5]and|6]were calculated by these methods. 
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